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PREFACE 


Meteorological  studies  suggest  that  technologically  feasible  oper¬ 
ations  might  trigger  substantial  changes  in  the  climate  over  broad  re¬ 
gions  of  the  globe.  Depending  on  their  character,  location,  and  scale, 
these  changes  raight  be  both  deleterious  and  irreversible.  If  a  foreign 
power  were  to  bring  about  such  perturbations  either  overtly  or  covertly, 
either  maliciously  or  heedlessly,  the  results  might  be  seriously  detri¬ 
mental  to  the  security  and  welfare  of  this  country.  So  that  the  United 
States  may  react  rationally  and  effectively  to  any  such  actions,  it  is 
essential  that  we  have  the  capability  to:  (1)  evaluate  all  consequen¬ 
ces  of  a  variety  of  possible  actions  that  might  modify  the  climate, 

(2)  detect  trends  in  the  global  circulation  that  presage  changes  in 
the  climate,  either  natural  or  artificial,  and  (3)  determine,  if  pos¬ 
sible,  means  to  counter  potentially  deleterious  climatic  changes.  Our 
possession  of  this  capability  would  make  incautious  experimentation 
unnecessary,  and  would  tend  to  deter  malicious  manipulation.  To  this 
end,  the  Advanced  Research  Projects  Agency  initiated  a  study  of  the 
dynamics  of  climate  to  evaluate  the  effect  on  climate  of  environmental 
perturbations.  The  present  Memorandum  is  a  technical  contribution  to 
this  larger  study. 

Rand's  position  ou  climate  and  weather  modification  studies  was 
asserted  in  its  publications  RM-3205-NSF  and  RM-5835-NSF.  The  approach 
to  understanding  the  consequences  of  climate  change  must  consist  of 
many  converging  paths;  and  for  reasons  of  safety  and  economy  this  ap¬ 
proach  must  take  advantage  of  the  growing  versatility  of  computers. 

The  present  Memorandum  is  one  of  a  series  of  numerical  investi¬ 
gations  of  the  wind-driven  oceanic  circulation,  and  is  intended  to 
broaden  our  knowledge  of  the  ocean/atmosphere  interactions  that  are 
important  to  weather  and  climate.  Other  Rand  publications  related 
to  the  present  Memorandum  include  RM-6110-RC  and  RM-6211-ARPA. 
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ABS TRACT 

From  the  time-dependent  numerical  solutions  of  a  wind-driven  homo¬ 
geneous  ocean  model,  a  negative  lateral  eddy  viscosity  of  the  order  of 
6  2  “1 

10  cm  sec  is  estimated  in  the  basin  interior,  where  the  flow  is 

dominated  by  transient  Rossby  waves.  These  waves  produce  a  systematic 

meridional  convergence  of  eddy  momentum  into  the  latitudes  of  mean 

eastward  current,  and  are  in  this  sense  analogous  to  the  Rossby  waves 

in  the  general  circulation  of  the  atmosphere.  This  regime  exists  in 

8  2  -1 

spite  of  the  formally  prescribed  lateral  viscosity  of  10  cm  sec  in 
the  model's  equations;  the  required  viscous  dissipation  (and  the  bulk 
of  the  meridional  momentum  transport)  apparently  occurs  in  the  standing 
waves  of  the  western  boundary  current  system.  These  results  suggest 
that  a  strong  gradient  (and  change  of  sign)  of  the  effective  eddy  vis¬ 
cosity  may  be  characteristic  of  an  ocean  with  meridional  boundaries. 
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I .  INTRODUCTI ON 


The  wind-driven  ocean  circulation  may  be  viewed  as  the  outcome  of 
a  hydrodynamical  balance  between  the  input  of  kinetic  energy  by  the 
surface  wind  stress  and  the  dissipation  of  energy  by  friction,  with  the 
inertial,  Coriolis,  and  pressure  forces  determining  the  local  charac¬ 
teristics  of  the  flow  for  a  particular  basin  geometry.  The  successful 
examination  of  the  steady  circulation  in  bounded  ocean  basins  thus 
dates  from  the  work  of  Stonmel  (1948),  who  introduced  a  bottom  friction 
proportional  to  the  (vertically  integrated)  transport,  and  Munk  (1950), 
who  introduced  a  lateral  eddy  viscosity  proportional  to  the  transport 
Laplacian.  Each  of  these  formulations  is  capable  of  simulating  a  west¬ 
ern  boundary  current  in  a  homogeneous  ocean,  and  one  or  the  other  has 
been  used  in  most  subsequent  studios  of  the  wind-driven  circulation. 
Although  the  mechanism  by  which  the  real  ocean  dissipates  the  energy 
of  the  large-scale  circulation  remains  uncertain,  it  appears  that  the 
assumption  of  a  lateral  viscosity  permits  the  more  realistic  descrip¬ 
tion  of  the  current’s  structure  near  the  western  shore,  in  that  the 
maximum  meridional  transport  occurs  close  to  (but  not  at)  the  boundary, 
and  a  realistic  lateral  countercurrent  is  permitted  (Stommel,  1965). 

A  continuing  criticism  of  the  (linear)  frictional  lormulatlon,  how¬ 
ever,  has  arisen  from  the  apparent  need  for  an  eddy  viscosity  A„  of  the 
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order  of  10  to  10  cm  sec  to  achieve  a  realistic  boundary-current 

ft 

width,  while  observational  studies  have  sometimes  yielded  eddy  viscos¬ 
ities  an  order  of  magnitude  smaller  (Stonmel,  1955;  Ichiye,  1957).  De¬ 
pending  upon  the  characteristic  eddy  scale,  viscosities  over  the  range 
10  to  10  cm  sec  have  been  reported  [see,  for  example,  Stommel  (1951, 
1955) ] ,  and  appear  to  be  in  accord  with  the  concept  of  a  power-law  de¬ 
pendence  of  the  horizontal  turbulent-diffusion  coefficient  upon  eddy 
size  (Stommel,  1949;  Ozmldov,  1968).  In  addition  to  this  scale  depen¬ 
dence,  the  present  note  suggests  that  the  circulation  may  be  characterized 

ft  —1/2  —1/3  1/3 

The  width  of  the  western  boundary  current  is  2tt3  8  A„ 

in  the  Munk  (1950)  theory  (where  8  is  the  northward  derivative  of  the 
Coriolis  parameter) ,  whereas  the  total  meridional  transport  of  the  cur¬ 
rent  is  independent  of  Ay. 
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by  a  strong  spatial  dependence  of  the  effective  eddy  viscosity,  taking 
on  negative  values  in  those  regions  dominated  by  transient  Rossby  waves. 


II.  ESTIMATION  OF  THE  EDDY  VISCOS ITY 


The  fluctuations  of  eddy  momentum  in  a  turbulent  flow  may  contri¬ 
bute  an  effective  (or  Reynolds)  stress  to  the  maintenance  of  the  time- 
averaged  flow,  and  may  therefore  be  regarded  as  giving  rise  to  an  eddy 
viscosity.  Considering  the  northward  eddy  transport  of  zonal  (eddy) 
momentum  as  a  stress  associated  with  the  mean  zonal  current  u,  we  may 
thus  write 


u  v 


-A 

3y 


(1) 


where  u'  and  v'  are  rvstward  and  northward  velocity  departures  from  the 
time-average  denoted  by  (  ) ,  and  y  is  the  northward  (cross-stream)  co¬ 
ordinate.  The  net  zonal  viscous  force  due  to  this  stress  is  thus 


f7<uV)  <2) 

-  —  2 

assuming  (u'v1)  and  (u')  to  be  independent  of  the  eastward  coordinate 
(and  depth),  and  assuming  zero  vertical  velocity  throughout. 

In  terms  of  the  eddy  viscosity  A^(y),  Eq.  (2)  leads  to 


dAH  3u  3ju 

dy  3y  ^  3y2 


-f7(uV) 


whose  solution  apart  from  an  arbitrary  constant,  is 


(3) 


(4) 


as  might  have  been  inferred  directly  from  Eq.  (1)  for  nonzero  3u/3y. 
In  the  case  of  vanishing  shear  3u/3y  *»  0,  Eq.  (3)  leads  to 
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These  solutions  thus  permit  the  estimation  of  in  terms  of  the  eddy 
momentum  flux,  with  the  relation  (5)  to  be  used  at  (or  near)  the  points 
of  maximum  or  minimum  mean  zonal  flow.  Here  we  have  in  effect  assumed 
that  u  (as  well  as  the  fluctuations)  is  uniform  in  the  zonal  direction, 
and  that  v  -  0.  A  more  formal  derivation  of  Eq.  (4)  may  be  found,  for 
example,  in  Proudman  (1953)  or  in  Neumann  and  Pierson  (1S66) , 


III.  THE  OCEANIC  MODEL  AND  ITS  SOLUTIONS 


The  evolution  of  the  horizontal  circulation  in  a  wind-driven  homo¬ 
geneous  ocean  model  of  uniform  depth  h  may  be  described  by  the  equations 


[~  (U2)  +  (UV  cos  <>)  -  UV  sin  ♦] 


9t  ah  cos  $ 


— 8— -  |f  +  fv  +  av2u  +  p"1^ 

a  cos  3X  s  o  X 


DV  S 

3t  f  h  cos  <t> 


[~  (UV)  +  |-  (V2  cos  4)  +  U2  sin  *] 


-  |f  -  fU  +  AV2V 
a  dtp  s 


H  +  JL_  + 

3t  a  cos  <J>  L  3X 


where  (U,V)  *  /,  (u,v)dz  are  the  vertically  Integrated  currents,  £  is 
-n 

the  free-surface  elevation,  <J>  and  X  are  the  latitude  and  longitude  on 

the  earth  of  radius  a,  f  is  the  Coriolis  parameter,  g  is  gravity,  A  the 

w 

horizontal  eddy  viscosity,  p  the  (uniform)  density,  t.  the  (zonal)  sur- 
face  wind  stress,  and  V*  the  spherical  surface  Laplaclan.  Here  bottom 
friction  has  been  neglected  and  S  is  a  parameter  which  determines  the 
degree  of  inertial  nonlinearity  through  the  assumed  vertical  profile  of 
the  horizontal  flow:  S  -  hJ|v|2dr/cT|v|dz)2.  Selected  parameters  in  a 

60°-longitude  basin  are  h  ■  400  m,  ■  -T  cos  ( 3n (<p  -  <j>  )/2(<i  -  *  )) 

-2  X  °8  2m 

with  T  ■  2  dynes  cm  ,  $  ■  18°  N,  4>  ■  66°  N,  and  A  ■  10  cm  sec  . 

o  m 

The  value  S  -  4  has  been  selected  as  representative  of  the  observed  ve¬ 
locity  profiles  in  the  Gulf  Stream  as  reported  by  Fuglister  (1963) ,  and 
introduces  only  a  modest  amount  of  nonlinearity.  The  solutions  were  de~ 
termined  from  a  centered  difference  analogue  of  Eqs.  (6)  to  (8)  from  an 
initial  state  of  rest;  further  description  of  the  model  and  details  of 
the  method  of  numerical  integration  are  given  elsewhere  (Gates,  1970). 
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Figures  1  and  2  show  the  30-to-60-day  n.aan  distribution  of  the 

vertically  averaged  currents  u  =  U/h,  and  v  ■  V/h,  respectively,  where 

(  )  denotes  the  time  mean  as  before.  Characteristic  of  these  solutions 

is  the  existence  of  a  near-equiiibrium  wester.i  boundary  current  system, 

together  with  a  series  of  transient  waves  wnich  dominate  the  interior 

circulation.  Also  shown  are  the  distributions  of  the  root-mean-square 

2  1/2 

(rms)  variabilities  u,v  ■  [(u',v')  ]  ,  where  (  )'  denotes  the  de- 

rms 

parture  from  the  time  mean  over  the  interval  30  to  60  days.  The  strong 
boundary  currents  in  the  west  are  relatively  steady  after  about  25  days, 
and  correspond  closely  to  those  jiven  by  the  linear  viscous  theory  of 
Munk  (1950) .  These  distributions  are  also  similar  to  those  shown  else¬ 
where  (Gates,  1970)  for  slightly  less  nonlinear  cases. 

In  the  interior  the  current  variability  is  associated  with  the 
westward  passage  of  successive  oceanic  Rossby  waves.  Except  for  a  slow 
amplitude  decay,  the  behavior  of  these  waves  away  from  the  western 
boundary  region  corresponds  reasonably  well  to  the  inviscid  theory  of 
transient  (free)  wave  motion  in  an  enclosed  basin,  with  the  reflection 
of  the  wave  energy  flux  from  the  western  shore  accounting  for  the  maxi¬ 
mum  variability  in  the  west-central  portion  of  the  basin  (Gates,  1969). 
In  spite  of  the  prescription  of  an  eddy  viscosity  in  the  basin-wide 
governing  Eqs.  (6)  and  (7),  the  interior  transient  solutions  are  there¬ 
fore  nearly  independent  of  the  prescribed  value  of  A. 
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IV.  THE  TRANSIENT  FLUX  AND  EDDY  VISCOSITY 

From  the  solutions  discussed  above,  the  mean  eddy-momentum  trans¬ 
port  u 1  v '  was  computed  and  was  then  averaged  zonally  across  the  basin 
(an  operation  denoted  by  []),  and  is  presented  in  Table  1.  In  view  of 
the  small  variability  in  the  western  boundary  currents,  the  value  of 
[u’v'j  is  principally  determined  by  the  interior  transients.  The  pro¬ 
file  of  the  mean  zonal  current  [u]  and  its  meridional  shear  3[u]/a3<|> 
given  in  Table  1  are  also  representative  of  the  interior  regime,  and  we 
may  therefore  use  these  quantities  to  determine  an  effective  interior 
eddy  viscosity  induced  by  the  transient  Rossby  waves.  On  the  assumption 
that  3u/3y  is  uniform  across  the  basin  interior,  the  eddy  viscosity  (A) 
may  thus  be  estimated  in  the  zonal  average  as 


Near  the  latitude  of  zero  mean  shear  3[u]/a3(|>  ■  0  (approximately  48°  N) , 
Eq.  (9)  was  replaced  by  Eq.  (5)  written  in  terms  of  zonally  averaged 
quantities  (see  Table  1) .  The  meridional  derivatives  in  these  calcula¬ 
tions  have  been  determined  from  data  (not  shown)  at  2°  north  and  south 
of  the  latitudes  indicated  in  the  table. 

The  meridional  profile  of  the  mean  zonal  current  [u]  shows  a  clearly- 
defined  eastward  flow  north  of  about  34°  N  in  response  to  the  eastward 
zonal  wind  stress  (which  reaches  a  maximum  at  50°  N) ,  except  for  the  mean 
westward  flow  north  of  58°  N,  which  is  a  recirculation  induced  by  the 
northern  boundary.  This  profile  gives  rise  to  a  nearly  uniform  positive 
lateral  shear  3[u]/ad$  over  most  of  the  southern  half  of  the  basin.  The 
profile  of  the  meridional  eddy-momentum  transport  [u'v' ]  also  shows  a 
clearly  defined  meridional  structure  that  is  quite  similar  to  that  of 
the  lateral  shear;  the  maximum  northward  flux  [u'v' ]  occurs  near  42°  N, 
close  to  the  maximum  positive  shear  at  44°  N.  The  profile  of  the  eddy- 
momentum  transport  thus  shows  a  broad  zone  of  flux  convergence  between 
42°  N  and  56°  N,  just  in  the  region  of  maximum  mean  zonal  current.  These 
data  suggest  that  the  eddy-momentum  flux  contributes  to  the  maintenance 
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ol  the  mean  zonal  circulation,  much  in  the  manner  that  the  mean  zonal 
atmospheric  flow  is  maintained  by  the  large-scale  migratory  atmospheric 
disturbances,  as  first  shown  by  Starr  (1951).  The  contribution  of  the 
transient  oceanic  Rossby  waves  to  the  zonal  momentum  balance,  however, 
is  small  compared  to  that  of  the  standing  zonal  waves  (i.e.,  the  west¬ 
ern  boundary  currents) ,  and  the  bulk  of  the  imposed  zonal  wind  stress 
is  balanced  by  a  pressure  torque  exerted  against  the  eastern  and  west¬ 
ern  boundaries  (Gates,  1970). 

The  effective  eddy  viscosity  A^  produced  by  the  transient  waves 

is  also  shown  in  Table  1,  as  found  from  Eq .  (9)  for  the  tonally  averaged 

statistics.  Over  the  latitudes  of  significant  northward  eddy  transport, 

6  2  —1 

the  eddy  viscosity  is  of  the  order  of  10  cm  sec  and  negative ,  in¬ 
dicating  that  the  transient  waves  serve  to  increase  (rather  than  de¬ 
crease)  the  mean  zonal  momentum.  In  terms  of  the  model’s  kinetic  energy 
balance,  we  may  also  view  these  data  as  representing  a  systematic  energy 
transfer  from  the  eddies  to  the  mean  flow.  Similar  features  are  also 
found  in  other  solutions  of  the  present  wind-driven  model  reported  else¬ 
where  (Gates,  1970).  It  is  interesting  to  note  that  the  observations 

of  Webster  (1965)  also  imply  a  negative  eddy  viscosity  of  the  order 
6  2  —I 

10  cm  sec  for  the  fluctuations  on  the  shoreward  side  of  the  Gulf 
Stream.  A  negative  eddy  viscosity  may  therefore  characterize  those  re¬ 
gions  in  the  ocean  which  are  dominated  by  transient  waves,  which  for 
the  present  model  includes  most  of  the  basin  interior.  Assuming  that 
the  necessary  lateral  viscous  dissipation  actually  occurs  in  the  bulk 
of  the  main  boundary  currents  (and  their  larger-scale  meanders),  a 
marked  spatial  variation  in  the  magnitude  of  the  eddy  viscosity  deduced 
from  velocity  fluctuations  may  thus  be  a  normal  feature  of  the  oceanic 
circulation. 
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V. _ CONCLUSIONS  AND  SPECULATION 


In  spite  of  the  prescription  of  an  overall  lateral  eddy  viscosity 
8  2  “1 

(A  =  10  cm  sec  )  in  the  wind-driven  model,  the  nearly  inviscid  tran¬ 
sient  solutions  in  the  interior  of  the  basin  are  responsible  for  a 

6  2  — 1 

negative  eddy  viscosity  of  the  order  10  cm  sec  .  The  oceanic  Rossby 
wave  regime  in  the  basin's  interior  is  tl.us  another  example  of  a  nega¬ 
tive-viscosity  phenomenon  (Starr,  1968),  and  is  in  this  respect  analo¬ 
gous  to  the  atmospheric  circulation.  The  oceanic  circulation  may  in 
fact  be  characterized  by  a  wide  range  of  effective  eddy  viscosity,  and 
the  reported  estimates  of  the  dissipation  on  selected  scales  and  in 
certain  regions  [such  as  those  of  Webster  (1965)]  do  not  necessarily 
invalidate  the  concept  of  a  frictionally  controlled,  wind-driven  cir¬ 
culation. 

In  the  bounded  basin  considered  here,  the  contribution  of  the  tran¬ 
sient  waves  to  the  momentum  transport  and  kinetic  energy  is  overshadowed 
by  that  of  the  standing  waves  of  the  western  boundary  current.  In  an 
ocean  without  meridional  boundaries,  and  hence  presumably  without  west¬ 
ward  intensification  (ignoring  topographic  effects),  the  transient  waves 
may  play  a  much  more  important  role  in  the  overall  circulation.  This 
may  be  the  case  in  the  circumpolar  seas  surrounding  Antarctica,  where 
it  seems  likely  that  the  persistent  wind  systems  would  excite  the  Rossby 
wave  modes.  Determination  of  the  role  of  the  oceanic  Rossby  waves  in 
this  and  in  other  oceans  of  the  world,  however,  awaits  the  availability 
of  suitably  long-term  observations  of  the  large-scale  flow. 
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